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Starting from a continuum description, we study the non-equilibrium roughening of a thermal re- 
emission model for etching in one and two spatial dimensions. Using standard analytical techniques, 
we map our problem to a generalized version of an earlier non-local KPZ (Kardar-Parisi-Zhang) 
model. In 2+1 dimensions, the values of the roughness and the dynamic exponents calculated from 
our theory go like a ~ z ~ 1 and in 1+1 dimensions, the exponents resemble the KPZ values for low 
vapor pressure, supporting experimental results. Interestingly, Galilean invariance is maintained 
althrough. 
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The subject of kinetic roughening and non-equilibrium 
growths, have been in the center of interest of far-from- 
cquilibrium physics for more than two decades now. This 
is mainly due to two reasons: on the one hand, due to an 
ongoing revolution in the world of micro-physics in recent 
years, the demand of the age is to understand and imple- 
ment the underlying mechanism associated [1]. On the 
other hand, they seem to correlate fields even as diverse 
as ecological growths, propagation of a crack-front, stock- 
market predictions, etc. [2]. Although the processes 
which have been probed so far, have mostly been con- 
cerned only with local effects, such as molecular-beam- 
cpitaxy (MBE) growth, conventional diffusive growths, 
etc, the importance of the non-local effects, have been 
known as early as the 1950's [3]. Later on, with the ad- 
vent of more sophisticated experimental techniques, non- 
linear effects involving physical vapor deposition (PVC) 
[1,4,5,6], sputtering techniques and associated growth 
and etching of plasma fonts have assumed a position of 
paramount importance. Whereas in standard MBE type 
of growths, the vapor atoms are targetted in a direction 
normal to the substrate, so that growth is decided by 
the local environment only, in case of shadowing growths 
by sputter deposition, vapor atoms are incident at ran- 
dom angles to the surface, so that non-local factors gain 
prominence in this case [7-11]. There have been several 
experimental follow-ups too of this sputtering mechanism 
[12-14]. 

The concept of shadowing effect in a sputtering growth 
(or etching) essentially arrived with the observation that 
thin films often exhibit " an extended network of grooves 
and voids in their interiors" [11] giving rise to columnar 
structures. The basic idea is the following. Since, in a 
sputtering growth (etching), particles are allowed to be 
deposited (deroded) on the surface from all possible an- 
gles at random, the rate of growth is taken to be propor- 
tional to the exposure angle 0(x), which is a function of 
the position of incidence of the incoming particle. Now, 
as the hills have greater exposure area, they receive more 
atoms than the valleys. Thus the hills continue to grow 



steeper compared to the depleted valleys, which natu- 
rally gives rise to an instability in the system. The idea 
has been very ingeniously, but intelligently related to the 
growth of the relatively larger stalks, in a grassy lawn, 
which suppress the growth of the shorter ones [11] and 
in the process giving rise to a rough contour. 

In the theoretical front, this phenomenon of shadow- 
ing growth (decay), or its partner, the thermal reemis- 
sion instability has inspired a series of works in 1+1 di- 
mensions [7,8,11,15-17] and in 2+1 dimensions [14,18,19]. 
The theoretical forays in fact started with the paper by 
Karunasiri, et al [7] where from a direct numerical in- 
tegration of the dynamical equation, they were able to 
show that the self-similarity of the contour, evident at 
small values of the diffusion constant, is modified by the 
growth of fiat films, beyond a critical height, as the value 
of the diffusion constant is increased. Taking clues from 
their arguments, Roland and Guo [15] went on to calcu- 
late the value of the roughness constant, in 1+1 dimen- 
sions (albeit in the context of a shadowing model) and 
further predicted that in the low temperature phase, the 
system resembles a KPZ universality class (in agreement 
with Karunasiri, et al [7]). This concept of non-local, 
shadowing effect was later modified in [9,11], where a 
net non-local flux was observed to give rise to the in- 
herent columnar structures found in experiments. Later 
on, the domain of 2+1 dimension was also probed with 
the advent of advanced numerical integration algorithms 
and Monte-Carlo simulations [18,19]. However, all these 
attempts, both in 1+1 and 2+1 dimensions, being pre- 
dominantly numerical, either through direct numerical 
integration of a fundamental Langevin-type equation, or 
through Monte-Carlo simulation, and all the more, giv- 
ing contradictory values of the exponents obtained by 
different groups, we ventured an analytical derivation to 
have a final say regarding the universality class of these 
type of sputtered mechanisms. In the process, we will 
see that our findings correlate the available experimental 
and numerical observations (of one of these groups) in 
2+1 dimensions and predicts scaling in 1+1 dimensions 
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too. 



With the assumption that the shadowing effect pro- 
vides the dominant instability in the system, we apply 
the non-local model proposed by Zhao, et al [14,18,19]. 
The model is given by 



Of 



dhij ' f) = uV 2 h(r, t) ± Jl + (Whf R(r, t) + V {r, t) 



and 



(V{f,t)r]{r',t')) = 2D 6 2 (F-r') 5{t-t') 



(1) 



(2) 



where the first term on the right hand side of eqn.(l) 
provides the diffusive relaxing mechanism for the grow- 
ing (or etching) surface and the last term signifies the 
collective effect of randomness in the system, taken to be 
a Gaussian noise. The middle term is the non-local, non- 
linear term detailing the effects of thermal reemission and 
is given by 



R{r,t) = s Q F Q (f,t) + s l F l {r,t) 



(3) 



Here sq is the zeroth order sticking coefficient and Si 
is generated due to the reemission mechanism [14]. Here 
we consider first-order thermal reemission, that is neglect 
the effects of Si (i > 1). Plugging again from the same 
reference and applying the same logic, we consider the 
flux of the m-th order particle at position r as F m (r,t) 
which is given by 



(4) 



F m+1 {r, t) = (1 - s m ) / Z(r, 7, t) F m {P , t) x 



ifippi .h) P(hp?, n') 
(r-r 1 ) 2 + (h- h'f 



dA' 



For our case of first-order reemission, we are concerned 
with m=0 and 1. Here h is the unit normal to the surface 
at r, pointing outwards, n' is the unit normal at r 1 and 
the unit vector connecting r and f 1 (see Fig.l). 
P(np/p, h') is the probability, per unit solid angle that 
the reemitted particle flies off along npip and is expressed 
as ^ Flr "'" [18]. Z(r,r' ,t) is equal to unity except when 
there is no line of sight between the surface elements 
at r and r and zero otherwise. The nonlinear factor 



1 + (Vft.) which is multiplied with i?(r, t), signifies the 
lateral growth (or etching, as the case may be) associated 
and the '+' and '-' signs as its prefix, refer to growth or 
etching respectively. In the following analysis, we will 
consider parameter values as in [19] (that is we will be 
dealing with etching due to sputtering). Thus, for our 
case, F = 4, s « and Si « 1. Also P(fip,f,n') = 
-SLsUl ; assuming thermally re-emitted flux, although this 
is more of a simplification [20] than exact truth. With the 
above description of the complete equation, we proceed 



to determine the dynamics in the 2+1 dimensional case. 
Later on, we will also discuss our results with reference 
to 1+1 dimensions, as well. 
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FIG. 1. Relative orientations of the unit normals at f and 
r* and the co-ordinate system described by them. 

Combining eqns.(l), (3) and (4) and taking ip as the 
angle between r and f 1 (see Fig.l), the dynamical etching 
equation reduces to 



^ « vV 2 h [1 + ^(V^) 2 ] Ftft) + V (f,t) 



(5) 



where 

Fi(r,t) 



4cos6» sm[0 + 6'] 



(r-r 1 ) 2 + {h-h'f 



1 + (W(F,0')) r'dr'dO' 



(G) 



where 8 = angle between hpf and n' — cf> + -0 as in 
Fig.l and 9' is again defined as in Fig.l. In arriving 
at eqns.(5) and (6), we have deliberately chosen n as 
one of the axes in the two dimensional plane, to simplify 
calculations. This can be done, since on the average this 
holds true. Also the standard lateral growth assumption, 
\Vh\ < 1 has been employed too. This Fi(r,t) can be 
further reduced to 



i< COS0 > 



dr 



, y -r\[i + \{d r/ h'f 



[h 



(7) 



where L is the size of the system. It is important to men- 
tion here that in deriving eqn.(7) from eqn.(6), we have 
used the mean-value theorem, since ir/2—5 < 6' < ir/2+S 
(d is an angular strip around h), the range being evi- 
dent from Fig.l. The "«" sign justifies the fact that 
we have taken a mean-valued average, represented by 
"<>" around the h-axis, thereby removing < cos# > out- 
side the integral as a first-order approximation. Simpli- 
fying further, we arrive at the analytically tractable form 
of F\(r,t), as given below: 
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8< cosfl > 2 



L 



Ud r ,h'f 



(8) 



In arriving at the above equations, we have put on 
a very standard assumption for any non-local model 
that the height difference (h — h'), calculated between 
any two points r and r* of the growing surface should 
be much smaller than their distance of separation, i.e. 



\h-h'\ < 



| , a basic property expected of any non- 



local process. 

With this assumption and the mean-valued average 
done beforehand, the equation of motion now becomes 



dh 8<cos0> 5 
at 7r 



4< cos# > l 



dr'- 



(d r ,h(r')f + n(r 7 t) (9) 



Now, we try to look at the possible large time, long 
distance behavior of the system. We can easily see that 
the KPZ part [21], constituting the second term on the 
R.H.S. of the above equation will vanish as the system 
size is taken to be sufficiently large. In deriving the above 

-» 2 

form, terms higher than (S7h) order have been neglected. 
The final equation now looks like 

— ^vV 2 h+ dr'cj)(r,r')\dr'h\ 2 +rj(r,t) (10) 



where 



(r, r') 



4Att 



(11) 



A = Ao < cos 8 > is an adjustable coupling parameter, 
such that we will later put Ao equal to unity. The fact 
that the assumptions employed above are perfectly trust- 
worthy, can be cross-checked from the fact that eqn.(ll) 
maintains translational invariance which was an impor- 
tant feature of our starting eqn.(4). 

Eqn.(10) can be easily mapped to the phenomcnolog- 
ical equation considered in [22]. The only trick lies in a 
suitable wave- vector representation of the effective long- 
range potential 4>(r,r') in our case. Obviously, this can- 
not be a simple plug-in from the earlier equation of mo- 
tion [22], since, here, the interacting potential is appar- 
ently a multi-valued function. To progress further, we 
move on to the wave-vector representation of this inter- 
acting potential which is given by the scaled relation 



Here the scaling function looks like 



dX Xe 



-iX 



Ye 



-iY 



(Y-&X) 



(12) 



(13) 



Considering the scaling ansatz 



we get 



cj>(k,k') w A 



(14) 



(15) 



and our job now is to evaluate the definite scaling be- 
havior for f(k, k') by the evaluation of a number for r\ 
from eqn.(13) [24]. Applying simple Laplace transform 
and going through the standard steps, it is easy to see 
that the dominating contribution of the double integral 
in eqn. (13) implies that rj = 1 [23] and this gives the 
value 



4>(k,k') w A 



(16) 



i. e. the major contributing part of the potential is 
effectively reduced to a single variable mode. Now, we 
can simply plug-in results from ref. [22] and write down 
the dynamic exponent z as 



where 



K 



z = 2 + K 



-24/23 



1.04 



(17) 



(18) 



for our case [24]. One obvious point to be noted here is 
the fact that owing to the Galilean invariance of eqn. (9), 
we can easily see that 



z = 2 



(19) 



and interestingly enough, the general tendency of the 
system is to flow towards a short-ranged fixed point (the 
long-ranged fixed point comes out to be unphysical with 
the specific parameter values, for our particular case). 
This effect, as we will see, holds sway in 1+1 dimensions 
too, where the system flows towards the KPZ fixed point. 
Combining the last two equations, we get 



a 



-K 



(20) 



Thus the critical exponents come out as 
22 

z = — = 0.96 
23 

24 

a =23 = L ° 4 

/? = - » 1.08 

z 



(21) 



i. e. a s (3 rj z m 1 in reasonable agreement with 
experimental and numerical findings [14,18,19] (experi- 
mental values are: a = 0.96 ± 0.06, (3 = 0.91 ± 0.03 and 
z = 1.05±0. 08), within experimental error bars. The fact 
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that the theory (and also experiment [18]) predicts a w 1 
indicates that the effects of overhangs might be marginal 
(pg. 110, ref. [1]). Also to be noted is the invariance of 
the Galilean identity a + z = 2. Before concluding this 
portion, it must be mentioned that for the opposite sce- 
nario, i. e. growth under first-order thermal reemission, 
an identical analysis as above shows immediately that 
now the reduced dynamical equation has a form nearly 
the same as that in eqn.(10) but with a negative non- 
local potential. This automatically suggests that due to 
the attractive nature of this potential, the growth finally 
stops at sufficiently large times ("smoothens") and j3 rs 
[18]. Interestingly, we find that even without thermal 
reemission, this marked change in the scaling properties, 
depending on whether it is a growth or an etching process 
has been discussed elsewhere [25] also. 

For the 1+1 dimensional case, we follow exactly sim- 
ilar lines, the only modification being the consideration 
of 6' = and 9 = o or n (depending on growth or de- 
cay, respectively) in eqn.(6). Thereafter, proceeding like- 
wise, the dominating long-ranged part comes out to be 
v(r) dr'(d r >h) 2 , with v(r) w |s Thus in the large time 
limit, as r — » L, we see that the system approaches the 
conventional KPZ fixed point and naturally the expo- 
nents too resemble the KPZ universality class, which can 
be looked upon as sort of an analogy with the shadow- 
ing case [7] . To avoid unnecessary repetition of identical 
calculations, as in the 2+1 dimensional case, we have 
neglected any further details in 1+1 dimensions. 

All said and done, however, there is still one open 
question which needs to be resolved. This is the fact 
that inspite of both the available short-ranged and long- 
ranged fixed points in the 2+1 dimensional case, the sys- 
tem chooses the short-ranged fixed point (an alternative 
statement that there is Galilean invariance in the system, 
since the other fixed point basically gives an unphysical 
picture with a < 0) although the shadowing effect fun- 
damentally remains a non-local contribution. This seems 
to suggest that whenever we are talking about non-local 
interactions, it does not necessarily mean that the long- 
ranged structure should control the associated dynamics. 
Instead the short-ranged part of the contribution might 
also take the upper hand, though, obviously depending 
on the type of interaction we are considering. The issue 
seems to demand further studies. As an adjoindcr, we 
would like to mention that the 1+1 dimensional situa- 
tion, being basically dominated by the KPZ fixed point, 
no such complexity arises over there. 
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